On Maxwell fluid with relaxation time and viscosity depending on the pressure 
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i — We study a variant of the well known Maxwell model for viscoelastic fluids, namely we consider the Maxwell 

^ffuid with viscosity and relaxation time depending on the pressure. Such a model is relevant for example 

ifjn modelling behaviour of some polymers and geomaterials. Although it is experimentally known that the 

^triaterial moduli of some viscoelastic fluids can depend on the pressure, most of the studies concerning the 

^motion of viscoelastic fluids do not take such effects into account despite their possible practical significance 

(^ui technological applications. Using a generalized Maxwell model with pressure dependent material moduli 

0\ve solve a simple boundary value problem and we demonstrate interesting non-classical features exhibited 

T"by the model. 
> 
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1. Introduction 

The Maxwell fluid model was originally developed by Maxwell [1] to describe the elastic and viscous 
response of air. Nowadays, it is however, frequently used to model the response of various viscoelastic 
fluids ranging from polymers — see for example Ferry [2] — to the Earth's mantle — see for example Cathles 
[3]. In the present paper, we study an important generalization of the original model due to Maxwell, 
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namely we consider a model with pressure dependent material moduli. While the model with a pressure 
dependent material moduli has important technological ramifications, little is known about the qualitative 
or quantitative features related to the model. In fact, there is no careful analytical study with regard to 
mathematical questions concerning existence and uniqueness of solutions for such fluids. Even within the 
context of solutions to initial-boundary value problems, there is no systematic study when all the material 
moduli are pressure dependent. Given the possible usefulness of this model, it is surprising that there are 
no such studies and the analysis carried out here addresses this lacuna. 

Bridgman [4] in his pioneering experiments in high pressure physics reported (for many organic fluids) a 
significant dependence of the viscosity on the pressure and articulated the need for the dependence of the 
material moduli on the pressure. Models of fluids with pressure dependent viscosity are nowadays frequently 
used to describe the behaviour of fluids in many applications, for example in lubrication theory, see Neale 
[5] and Gwynllyw et al. [6]. Extending to viscoelastic fluids, we can ask whether the material moduli 
of viscoelastic fluid models also exhibit dependence on pressure, and what effects can be described if we 
use pressure dependent material moduli. Viscoelastic material moduli depending on the pressure has been 
reported for polymers, see for example Singh and Nolle [7], McKinney and Belcher [8] and the literature 
stemming from these papers, as well as for geomaterials, see for example Weertman et al. [9], Ivins et al. [10] 
and Sahaphol and Miura [11]. 

The question on pressure dependent viscosity and/or relaxation time is especially interesting, for example, 
with respect to applications in geophysics, since the material of the Earth's mantle is subject to a wide range 
of pressures. If we consider a material stratified due to the influence of the gravitational force, then because of 
pressure dependent material moduli, we would be dealing with a body with material moduli depending on the 
vertical coordinate. A similar situation occurs in geophysical applications and usually it is assumed that the 
body (in this case the Earth's mantle) is composed of a number of layers of materials with constant material 
moduli — a paradigm introduced in papers by McConnel [12, 13]. The approach based on material moduli 
dependence on the pressure and consequently (in the case of a stratified material) on the vertical coordinate 
provides an alternative approach to the problem. In such an approach, the material moduli continuously vary 
with the depth, in contrast to McConnel [13]. Although geophysicists and polymer engineers are aware of the 
possibility of pressure dependent material moduli, this dependence is invariably ignored in studies concerning 
the dynamics of these materials. If we take into account that geologists often try to refine their models for 
the viscosity by considering advanced non-Newtonian models, for example a power-law type viscosity — see 
for example Weertman et al. [9], Wu and Wang [14] — one wonders why important physical phenomenon such 
as material moduli depending on pressure is not considered. 

In the present paper we would like to explore potential benefits of models based on the assumption 
that the material moduli are pressure dependent. We consider a simple boundary value problem for an 
incompressible Maxwell fluid with the relaxation time and/or viscosity depending on the pressure to illustrate 
the consequences of the material moduli depending upon the pressure. The constitutive model therefore has 
the following structure: 




T = -ttI + S, 



(1.1a) 
(1.1b) 
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where T denotes the Cauchy stress tensor, 7r is the pressure 2 . 

fi(ir) = HoO- + Pn), (1.2a) 

/i(7r) = noe?*, (1.2b) 

A(tt) = A (1 + 77 r) , (1.2c) 

A(tt) = A e^, (1.2d) 

where Ao, /J>o, P and 7 are given constants. Viscosity in the form (1.2a) and (1.2b) is used in studies dealing 
with viscous fluids with pressure dependent viscosity, see for example Hron et al. [17] and Srinivasan and 
Rajagopal [18]; models (1.2c) and (1.2d) are counterparts of (1.2a) and (1.2b) with respect to modulus A. 
From the theoretical point of view, all the models that are considered fit into the thermodynamical framework 
developed by Rajagopal and Srinivasa [19] (see also Rajagopal and Srinivasa [15]). 

The boundary value problem we are going to solve using the model (1.1) is a classical one; it is a 
variant of the second Stokes' problem introduced by Stokes [20]. The fluid is confined between two parallel 
plates, the bottom (top) plate is at rest and the top (bottom) plate is moving with time-periodic velocity 
V = V cos(wt)e £ . The pressure on the bottom plate is fixed to n and a specific body force — the gravitational 
force — is acting on the fluid, the gravitational acceleration is denoted as g. (See Figure 1.) We call the 
problem a variant of the second Stokes' problem, since Stokes [20] considered the problem in a half-space 
above an oscillating plate and without the presence of the gravitational force. A setting identical to ours was 
previously studied by Srinivasan and Rajagopal [18] in their paper concerning viscous fluids with pressure 
dependent viscosities. 

There is a further need for the careful assessment of the model by considering initial-boundary value 
problems that have technological relevance and also those that can be used to correlate the model with 
experiments. From the point of geophysical applications another simple problem that is worth studying 
using the model (1.1) is a variant of problem introduced by Haskell [21] (response of a fluid on the removal 
of load), and isostatic adjustment in incompressible, nonrotating and self-gravitating spherical planet, a 
problem introduced by Love [22] — see also discussion in Wolf [23]. 



V cos (ut) e$ 



A e fi 



Figure 1: Problem geometry. 



2 In this paper, we shall refer to the Lagrange multiplier ir as pressure. For the model (1.1), the Lagrange multiplier is not 
the same as the mean normal stress, and if one defines the mechanical pressure as being the mean normal stress, then 7r is not 
the mechanical pressure. For further discussion see Rajagopal and Srinivasa [15], Rajagopal [16]. 
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2. Governing equations 

We now turn our attention to solving the boundary value problem. Equations governing the motion of 
the material are 



dv 

9 ~dt = dlv T + 
div v = 0, 



(2.1a) 
(2.1b) 



where v denotes the velocity, p is the density b = — ge^ is the specific body force and T is given by 
constitutive model (1.1). No-slip boundary conditions on the top and bottom plate read — for the oscillating 
top plate — v| = and v| h = V cos(uit)ert, and vice versa for the oscillating bottom plate 3 . Since we are 
dealing with a simple geometry, the system (2.1) can be greatly simplified under the assumption of parallel 
flow. 



2.1. Governing equations for parallel flow 

Let us assume that the velocity field and the pressure field have the form v = v x (y,t)ex, n = n(y) 
respectively, and let us suppose that 8 = S(y, t). It follows that [Vv] v = 0, [VS] v = and 



L = 



® vX 

dy 
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Q ® vX 
dy 

dv x n 
. dy U 



s = 



2 dv^ 



dt 

®^xy dv x O 

dt dy ^yy 



dS t< 



dv x q 

Thi^xy dt dy D yy 

dt 



The balance of linear momentum (2.1a) reduces to 

dv x dS^ 

1 dt 



= - 



dy 

dn dS„ 



+ 



dy dy 



- 99, 



and (1.1b) reduces to 



S ££ + \(n)-^-2\(n)—S, § =0, 
dS~ ~ dv* 



dt 



dy yy 



dv x 
dy 



(2.2a) 
(2.2b) 



Obviously, if we are interested in a time-periodic solution, we need to fix 
solve (2.2b) to get a formula for the pressure 

7T = -pgy + TTO- 



(2.2c) 
(2.2d) 
(2.2e) 

0. If = 0, we can then 

(2.3) 



3 For the sake of brevity we will, however, describe a solution procedure only for the top plate oscillating, and we will only 
report the results for the case when the bottom plate is oscillating. 
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Finally, the system reduces to 



p w = if- < 2 - 4a > 



dSz, fl , fiv x 



S xx + - 2A(vr)— = 0, (2.4b) 



S*q + = Mt)-^", (2.4c) 

where X(ir) and are, by virtue of (2.3), known functions of y. Equations (2.4a) and (2.4c) are linear 
equations for v x and 5-- and these equations can be solved independently on (2.4b). The system has to be 
completed by supplying boundary conditions, in our case (oscillating top plate) the boundary conditions for 
the velocity read 

v% =h = V cos (ut), v\ =Q = 0. (2.5) 

2.2. Dimensionless governing equations for parallel flow 

Let us develop a dimensionless version of the governing equations (2.4) before we proceed to solving the 
problem. Let h be the characteristic length and V the characteristic velocity, then system (2.4) reads 

dv x * 1 dS%- t x 

(2.6a) 



dt* Re dy* ' 

dt* dy 



S xx + WeA*(7r*) ( - ) = 0. ( 2 - 6b ) 



^ + WeA>*)^ = //K)^, (2.6c) 

where Re = is the Reynolds number (the ratio of inertial forces to viscous forces), We = ^f- is the 

Weissenberg number (the ratio of the relaxation time of the material and a characteristic time scale), and 

the star denotes dimensionless variables v x * = y-, S^- = , S^- = y* — |, t* — Functions A*(7r*) 
and £t*(7r*) are given by formulae 

H*(n*) = (1 + /?*tt*) , (2.7a) 

fi*(n*)=er* k , (2.7b) 

A*(7r*) = (l + 7 V*), (2.7c) 

A*(7r*)=e^\ (2.7d) 



where (3* = /37r , 7* = 77r and 

7T* = -Ify* + 1, (2.8) 

where IT = ^ is a dimensionless version of (2.3). Clearly, the dimensionless version of boundary condi- 
tions (2.5) reads 

v x * =0, v £ * =cos(wV). (2.9) 

Hereafter, we will use only dimensionless variables and we will therefore omit the star denoting the dimen- 
sionless variables. Furthermore, in what follows we will freely vary all the dimensionless parameters that are 
present in the equations. The parameter values and their mutual relationship do not necessarily correspond 
to a realistic setting; the main aim of the present study is to illustrate the trends that correspond to particular 
variations of the parameters. 
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3. Solution to the governing equations 

Due to the linearity of the governing equations and the nature of the boundary conditions (2.5), we seek 
a solution of the form 

v x (y,t) = v x (y)e i "\ S, y (y,t) = S, y (y)e^, (3.1) 

where v x (y) is a complex function, solution to (2.6) is then the real part of v x (y)e luJt and S £y (y)e luJt respectively. 
Substituting (3.1) to (2.6) leads to system 

iuv* = J-^-, (3.2a) 
Re dy ' K J 

(l + kjWe\(*))S j 0= t M(*)^, (3.2b) 

that can be reduced to a single ordinary differential equation for the velocity component v x , 

£ _ 1 d / jxfr) dv*\ 

Redy Vl + iwWeA(7r) dy J ' K ' J 

that has to be solved subject to boundary conditions 

H {^) \y=0 = °' M (^ e ^) U = C ° S (^)' ( 3 - 4 ) 

where 9ft denotes the real part of the corresponding expression. The boundary conditions must be satisfied 
for all t. Equation (3.3) with boundary conditions (3.4) can be rewritten as a system of ordinary differential 
equations for the real and imaginary part of v x . Let and v x m denote the real and imaginary part of v x , 
then (3.3) reads as follows 



-uRev x m = - 



uRev Re = - 



A«M dvL , M*>WeA(7r) dvf m 



+ 



1 + (coWeX(n)y dy 1 + (wWeA(7r)) 2 dy 
^(Tr)wWeA(Tr) dv| e //(tt) dvf m 



+ 



(3.5a) 
(3.5b) 



1 + (wWeA(yr)) 2 dy 1 + (wWeA(7r)) 2 dy 
Using (3.1), the required velocity is given by 

4e = 4c cos(wt) - vf m sin(wi), (3.6) 
and the boundary conditions (3.4) reduce to 

4c|, =0 = 0, v x Rc \ y=1 = l, vf m \ y=o = 0, vf m \ y=1 = 0. (3.7) 
In addition, the vorticity only has a component is the z-direction whose value is given by 

( rot vf = = -^cosM) + ^sin(wf). (3.8) 
dy dy u/y 

5. J. Analytical solution to the governing equations 

Let us now solve (3.3) subject to boundary conditions (3.7). Although we will not be able to get a 
closed analytical formulae for the solution, we can, in most cases, give analytical formulae wherein certain 
integration constants have to be found numerically. 
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3.1.1. Case \(ir) = Ao, ^(n) = fio 

In the classical case (3.3) reduces to 



* 1 1 d 2 v x . , 

lUV = ReTT^We^ (3 ' 9) 

and the solution reads 

v*(y) = C ie iiBy + C 2 e- iiBy , (3.10) 



where B = A/cuRe (1 + iu/We), and Ci and C 2 are arbitrary complex constants. In this case it is in principle 
possible to apply the boundary conditions (3.7) and obtain values of C\ and C 2 as combinations of elementary 
functions of parameters Re, We and uj. (A solution to the classical second Stokes' problem in the infinite 
domain was given by Fetecau and Fetecau [24].) Also, note that in this case the velocity field is not influenced 
by the gravitational force; however, the gravitational force affects the pressure field. In the cases that follow, 
the gravitational force affects both velocity and pressure fields. 

3.1.2. Case X(tt) = A (1 + 77r) ; (j,(tt) = fi 

If X(tt) is given by (2.7c), then using (2.8) we can substitute for the pressure and get the relaxation time 
as a function of y, X(y) = 1 + 7 — jlly . Denoting y — 1 + iuWeA(7r), and using y as a new independent 
variable, v x {y) = v x (y), we see that (3.3) can be rewritten as 

V — 5 rr + i ~y 2 v x = 0. (3.11) 

df dy w (We 7 n) 2 V ; 

This equation is of form x 2 y" + axy' + (bx n + c) y = 0, b 7^ 0, n 7^ 0, for which the solution (see Polyanin 

and Zaitsev [25]) has the form y = x 1 ^ {c^ u (j^Vbx^ + C 2 Y U (j^Vbx^^j, where v = ^a/(1 - a) 2 - 4c, 

C\ and C2 are arbitrary constants, and J u (z) and Y u (z) denote Bessel functions of the first and second kind. 
Solution to (3.11) therefore reads 

v%y) =y (cy f {^By^j +C 2 Y § Qi'Sy')) , (3.12) 

where B = ^J~^ , y = 1 + iwWe (1 + 7 — 7 Ily) , and C\ and C 2 are arbitrary complex constants that can 

be fixed by applying the boundary conditions (3.4). Since there is no general formula for splitting a Bessel 
function of a general complex argument into its imaginary and real parts (see for example Watson [26]), 
the system of linear algebraic equations arising from the application of the boundary conditions contains 
factors — the Bessel functions evaluated at y = and y = 1 — that must be found numerically. An algorithm 
for numerical computation of Bessel functions of the complex argument is discussed for example by Amos 
[27]. Formula (3.12) is worth of considering in the case when it is necessary to have a kind of analytical 
expression for the solution. If one however needs only a numerical solution, it is simpler to directly solve the 
system (3.5), see Section 3.2. We will face the same situation in the remaining cases. 

3.1.3. Case A(tt) = A 0; = // (1 + f3n) 

If /i(7r) is given by (2.7a), then using (2.8) we can substitute for the pressure and get the viscosity /i(7r) 
as a function of y, n{y) = 1 + (3 — (3Hy. Denoting y — (1 + /? — ^Uy)^, and using y as a new independent 
variable, v x {y) = v x (y), we see that (3.3) can be rewritten as 



^ 2 d 2 v x Av x 4iwRe(l + iu/We) ^ £ 



f-HT + 7777^2 -V v x = 0. (3.13) 
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This equation is again of form x 2 y" + axy 1 + (bx n + c) y = 0, b ^ 0, n ^ 0. Solution to (3.13) therefore reads 

vHv) = V (CiJo (i f 5y) + C 2 Y , (3.14) 

where B = -j^^uRe (1 + io>We) , y — (1 + /3 — /?Ih/)2 ? and Ci and C 2 are arbitrary complex constants that 
can be fixed by applying the boundary conditions (3.4). 

3.1.4- Case A(7r) = A e 77r ; = /x 

If A(7r) is given by (2.7d), then using (2.8) we can substitute for the pressure and get the relaxation time 
X(n) as a function of y, X(y) = e 7(1 ~ n?/ ). Denoting y = iwWee 7 ^™ 11 ^, and using y as a new independent 
variable, v x (y) = v x (y), we see that (3.3) can be rewritten as 



f(l+y)^r- f^r -iB(l + y) v x = 0, (3.15) 
dy dy 



where B = w ^^ yT ^ ■ Equation (3.15) to the best of our knowledge, does not have an analytical solution. It 
is however possible to find a solution in terms of power series, but it would bring an additional complication 
to those discussed in Section 3.1.2, and we will therefore not give a solution in terms of power series. 

3.1.5. Case A(tt) = A , n(ir) = /i e^ 

If fj,(n) is given by (2.7b), then using (2.8) we can substitute for the pressure and get the viscosity /x(7r) 
as a function of y, /x(7r) = . Denoting y = (/i(7r))~ 5 = e~^^~^ Uy \ and using y as a new independent 

variable, v x (y) = v x (y), we see that (3.3) can be rewritten as 

„ 2 dV Av x 4iwRe(l + iwWe)^ 4 n fo . a . 

V — o- - y—r~ o -y v =0. (3.16) 

df dy (fill) 2 1 J 

This equation is again of form x 2 y" + axy' + (bx n + c) y = 0, b ^ 0, n ^ 0. Solution to (3.16) therefore reads 

vHv) = y (CiJi (i f 5y) + C 2 Y l , (3.17) 

where B = -|j^/uRe (1 + iu/We), y = e~^~^ Uy \ and C 1 and C 2 are arbitrary complex constants that can 
be fixed by applying the boundary conditions (3.4). 



3.2. Numerical solution 



For the sake of convenience, we shall denote - , ,, a by A and t^t^ttt^ by £>, and the terms in 

' l+(yWeA(7r))^ l+(u>WeA(7r)) j ' 

the square brackets of the right hand side of (3.5a) and (3.5b) by J 7 and Q respectively. Then, (3.5) can be 
re-written as a system of four first-order ordinary differential equations as follows 



dy 

dg 

dy 

dvj e _ AT-BQ 
dy ~ A 2 + B 2 ' 

dvf m _ AG + BT 
dy ~ A 2 + B 2 ' 



-uRevf m , (3.18a) 
uRevi c , (3.18b) 

(3.18c) 
(3.18d) 
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The system of ordinary differential equations (3.18) along with the boundary conditions (3.7) was solved 
using the solver bvp4c in MATLAB, for various cases of /i(7r), X(n). 

Plots of the solutions for various parameter values are given in Figures 2-14. Parameter values were 
chosen to demonstrate trends induced by their change and do not necessarily correspond to parameter values 
encountered in real situations. Figures 2-11 display solution for oscillating top plate, while Figures 12-14 
display the solution when the bottom plate is oscillating. 

Let us now consider the case when the top plate is oscillating. Snapshots of velocity profiles for various 
models are shown in Figure 2. Models (1.2a), (1.2c) and (1.2d), for given parameter values, lead to very 
similar velocity profiles, see Figures 2a, 2c-2e. Model (1.2b), however, for the given parameter values, shows 
a substantial departure from the classical model (constant viscosity and constant relaxation time), compare 
Figure 2a and Figure 2b. The difference between the velocity profiles is clearly visible in Figure 3a where 
we compare, at a particular time instant, the classical model and model (1.2c) for various values of 7. It is 
obvious that model (1.2c) allows the oscillations induced by the top plate to propagate up to the bottom 
plate. This holds for the other non-classical models as well, see Figure 3b and Figure 4, but model (1.2b) 
leads to strongest oscillations near the bottom plate. Moreover, model (1.2b) leads to the creation of layers 
where the velocity rapidly changes its direction with respect to the vertical coordinate. 

The fact that (1.2b) leads to the propagation of oscillations up to the bottom plate is also well documented 
by plots of vorticity given in Figures 7, 8, 9 and Figure 11. (The vorticity has, in the case of parallel flow, only 
one nonzero component, rot v = — ^-e^.) Heavily oscillating layers are generated in the fluid (1.2b) especially 
for high Weissenberg number (We), Reynolds number (Re) and, of course, 7; see Figure 5 and Figure 6 for 
the comparison of velocity profiles, and Figure 10 and Figure 11 for the comparison of vorticities. Obviously, 
with sufficiently large We, Re and 7, the model leads to a motion where the vorticity is concentrated near 
the bottom plate, although it is the top plate that forces the material to move, see Figure 11c. 

Similar effects can be observed in the case where the bottom plate is oscillating and the top plate is at 
rest, see Figures 12-14. Model (1.2b) again leads to the creation of highly oscillating layers near the bottom 
plate. On the other hand, model (1.2d) with increasing j3, leads to nearly vanishing velocity gradient near 
the bottom plate and steep velocity gradient near the top plate, see Figure 14a; one can say that the material 
develops a boundary layer close to the boundary opposite to the boundary inducing the motion (compare 
cases (3 = and j3 = 8)! 

4. Conclusion 

We have studied the time dependent flow of a generalized Maxwell model that can be used to describe 
viscoelastic materials in which the material moduli depend on the pressure. Such models are relevant in 
modelling behaviour of polymers and geomaterials. Using the generalized model we have solved a simple 
boundary value problem — a variant of the well known Stokes' problem for the flow induced by an oscillating 
plate. The solutions were found numerically, and for some particular cases we have also derived analytical 
expressions for the solution. We have shown that the pressure dependent material moduli reduce to material 
moduli that are continuously dependent on the vertical coordinate, and that the non-constant material 
moduli, compared to the classical Maxwell model, have substantial impact on the dynamical behaviour of 
the fluid. The departures from the behaviour predicted by the classical model are the most distinctive for 
the model with constant viscosity and relaxation time depending exponentially on the pressure. If the latter 
model is used instead of the classical one for the case when the top plate is oscillating, then the oscillations 
induced by the moving top plate can propagate up to the bottom plate, and the vorticity can reach its 
maximum near the bottom plate; these effects however cannot be captured by the classical Maxwell model. 

Although the fact that material moduli for viscoelastic fluids can depend on the pressure is known, it 
is rarely used in dynamical considerations, even if one expects that the arising pressure to be significantly 
different at different parts of the body. The present paper provides a simple illustration of the effects that 
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-1 -0.5 0.5 1 



v x 

(d) /i exponentially depending on pres- 
sure with /3 = 6 and A constant 




-1 -0.5 0.5 1 



(e) jU linearly depending on pressure 
with (3 = 6 and A constant 



Figure 2: Velocity profiles for various cases of relaxation time (A) and viscosity (fi) depending on pressure 
at different times, Re = 100, We = 0.1, w = 1, II = I. 




12 -0.1 0.1 0.2 0.3 0.4 0.5 




(a) A exponential, /i constant 



(b) A linear, [i constant 



Figure 3: Velocity profiles at cut — | for various 7 values when the relaxation time (A) depends on pressure 
exponentially and linearly. Viscosity (fx) is kept constant in both cases with Re = 100, We = 0.1, oj — 1, 

n = 1. 
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-0.05 0.05 0.1 0.15 0.2 0.25 0.3 0.35 -0.05 0.05 0.1 0.15 0.2 0.25 0.3 0.35 

(a) A constant, /i exponential (b) A constant, /i linear 

Figure 4: Velocity profiles at ujt = | for various j3 values when the viscosity (/x) depends on pressure 
exponentially and linearly. Relaxation time (A) is kept constant in both cases with Re = 100, We = 0.1, 
w = 1, n = 1. 













■~^TTT5* We = 1 




• We = 5 
J — ■ — We = 10 

1 



-1 -0.5 0.5 1 1.5 



(a) A exponential, [i constant, 7 = 
classical case 




-1 -0.5 0.5 1 



(b) A exponential, fi constant, 7 = 3 




-10 -5 5 10 



(c) A exponential, /i constant, 7 = 6 



Figure 5: Velocity profiles at cot — ^ for various We values for the classical case as well as when the relaxation 
time (A) depends on pressure exponentially. Viscosity (//) is kept constant in all cases with Re = 100, co — 1, 
IT = 1. 




(a) A exponential, /1 constant, 7 = (b) A exponential, \x constant, 7 = 3 (c) A exponential, [i constant, 7 = 6 
classical case 



Figure 6: Velocity profiles at ut = | for various We values for the classical case as well as when the relaxation 
time (A) depends on pressure exponentially. Viscosity (/i) is kept constant in all cases with Re = 1, u — 1, 
IT = 1. 



11 



(a) classical solution 



(b) A exponentially depending on pres- (c) A linearly depending on pressure 
sure with 7 = 6 and /x constant with 7 = 6 and \i constant 





(d) ti exponentially depending on pres- (e) fj, linearly depending on pressure 
sure with (5 = 6 and A constant with (5 = 6 and A constant 

Figure 7: Snapshots of vorticity profiles for various cases of relaxation time (A) and viscosity (//) depending 
on pressure, Re = 100, We = 0.1, to — 1, II = 1. 




(a) A exponential, /1 constant 



10 20 30 40 




(b) A linear, [i constant 



Figure 8: Vorticity profiles at ut = | for various 7 values when the relaxation time (A) depends on pressure 
exponentially and linearly. Viscosity (/i) is kept constant in both cases with Re = 100, We = 0.1, u — 1, 

n = 1. 
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(a) A constant, \i exponential 



(b) A constant, /i linear 



Figure 9: Vorticity profiles at cot = | for various /3 values when the viscosity (/x) depends on pressure 
exponentially and linearly Relaxation time (A) is kept constant in both cases with Re = 100, We = 0.1, 

oj = i, n = i. 




(a) A exponential, /i constant, 7 = (b) A exponential, fj, constant, 7 = 3 (c) A exponential, [i constant, 7 = 6 
classical case 



Figure 10: Vorticity profiles at ut = | for various We values for the classical case as well as when the 
relaxation time (A) depends on pressure exponentially. Viscosity (/i) is kept constant in all cases with 
Re = 1, u = 1, n = 1. 




(a) A exponential, /1 constant, 7 = (b) A exponential, \x constant, 7 = 3 (c) A exponential, /1 constant, 7 = 6 
classical case 



Figure 11: Vorticity profiles at cot = | for various We values for the classical case as well as when the 
relaxation time (A) depends on pressure exponentially. Viscosity (/i) is kept constant in all cases with 
Re = 100, u = 1, IT = 1. 
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(a) classical solution (b) A exponentially depending on pres- (c) A linearly depending on pressure 

sure with 7 = 6 and \i constant with 7 = 6 and \i constant 





(d) fj, exponentially depending on pres- (e) fj, linearly depending on pressure 
sure with (3 = 6 and A constant with (3 = 6 and A constant 

Figure 12: Snapshots of velocity profiles for various cases of relaxation time (A) and viscosity (/i) depending 
on pressure, with bottom plate oscillating and for Re = 100, We = 0.1, co — 1, II = 1. 




-1 -0.5 0.5 1 -0.2 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 



(a) A exponential, /1 constant (b) A linear, /i constant 

Figure 13: Velocity profiles at ut = | for various 7 values when the relaxation time (A) depends on pressure 
exponentially and linearly, when the bottom plate is oscillating. Viscosity (/i) is kept constant in both cases 
with Re = 100, We = 0.1, 00 = 1, II = 1. 
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(a) A constant, /x exponential 



(b) A constant, fi linear 



Figure 14: Velocity profiles at ut — | for various (3 values when the viscosity (/j) depends on pressure 
exponentially and linearly, when the bottom plate is oscillating. Relaxation time (A) is kept constant in both 
cases with Re = 100, We = 0.1, u = 1, n = 1. 

can be expected if the fact that the material moduli depend on pressure is taken into account in studying 
the motion of viscoelastic fluids. 
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